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ABSTRACT 


The estimation of planetary magnetic fields from observations of the 
magnetic field gathered along a spacecraft flyby trajectory is examined with 
the aid of generalized inverse techniques, with application to the internal 
magnetic field of Jupiter. Model non-uniqueness resulting from the limited 
spatial extent of the observations and noise on the data is explored and 
quantitative estimates of the model parameter resolution are found. The 
presence of a substantial magnetic field of external origin due to the 
currents flowing in the Jovian magnetodisc is found to be an important source 
of error in estimates of the internal Jovian field, and new models explicitly 
incorporating these currents are proposed. New internal field models are 
derived using the vector heliun magnetometer observations and the high field 
fluxgate observations of Pioneer 11, and knowledge of the external current 
system gained from the Pioneer 10 and Voyagers 1 and 2 encounters. 

INTRODUCTION 

Four spacecraft have thus far encountered Jupiter and permitted in situ 
investigation of many aspects of the Jovian system, including detailed 
observations of the inner magnetosphere. These spacecraft have explored 
different regions of the inner magnetosphere by virtue of their unique 
trajectories. The near-equatorial approach of Pioneer 10 to within 2.8 Jovian 
radii permitted the first observations of what has become known as the 
magnetodisc (Smith et al., 1974; Van Allen et al., 1974), a thin annular disc 
of tenuous plasma and charged particles encircling the giant planet to 
distances approaching 100 Jovian radii. Large scale azimuthal currents 
flowing in this magnetodisc, subsequently observed by the Voyagers 1 and 2 
spacecraft (Ness et al . , 1979a; Ness et al., 1979b; Bridge et al., 1979) lead 
to a substantial magnetic field of ext^jrnal origin throughout the entire inner 
magnetosphere. Tne high- inclination retrograde approach of the Pioneer 11 
spacecraft yielded measurements of the inner magnetosphere at high latitudes 
spanning a wide range of Jovian longitude. In addition, the close approach of 

Jovian radii (Rj) made this trajectory the most favorable for the 
purposes of estimating the internal magnetic field of Jupiter. Thus, models 
of the internal field based on the Pioneer 10 measurements (Smith et al., 

1974) were quickly supplanted by models based on the Pioneer 11 observations 
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(Davis and Smith, 1976; Aouna and N«ss, 1976). The near equatorial approaches 
of the Voyager 1 and 2 spacecraft, at the close approaches of >^5 and 10 Rj, 
have permitted more detailed studies of the external azimuthal current system 
(Connerney et al., 1981) but have thus far idded little to the present 
understanding of the internal magnetic field of Jupiter (Ness et al., 1979b). 

As the nunber of available Jovian internal field models grows it becomes 
increasingly important to evaluate the model non-uniqueness inherent in the 
estimation of planetary magnetic fields from flyby observations. While some 
of the proposed racKitls represented an intermediate stage of data processing 
and analysis, in general differences between the proposed models reflect 
differences in the number and kind of coefficients used in the basic models, 
differences in the actual observations, data intervals chosen for analysis, 
and the choice of weights applied to the observations in the least squares 
minimization process, A recent summary of the magnetic field observations of 
Jupiter (Smith and Gulkis, 1979) lists six spherical harmonic models of the 
Jovian magnetic field; five based on the Pioneer 10 and 11 vector heliun 
magnetometer (VHM) data (Smith et al., 197^; Smith et al., 1976) and one 
derived from the Pioneer 11 fluxgate magnetometer (FGM) observations (Acuna 
and Ness, 1976). These models (with the exception of P10-11 combined models, 
e>g.. Smith and Gulkis, 1979) all represent a ’good’ fit along one trajectory 
in the sense of minimizing the weighted or unweighted model residuals, but 
lead to substantially different model fields in other regions of space. 

The model non-uniqueness we refer tc arises from the spatial limitations 
of the available observations and the noise on the data. In the usual 
spherical harmonic representation of planetary fields, there exist certain 
linear combinations of the model parameters which lead to a very small 
magnetic field along the spacecraft trajectory, but a large magnetic field 
elsewhere (e.g., at the planet's surface). That is, observations along a 
single flyby trajectory are insensitive to certain combinations of parameters. 
This basic lack of information, inherent in the available observations, leads 
to an unavoidable model non-uniqueness. 

This paper concerns the estimation of planetary magnetic fields from 
flyby observations, in particular from the Pioneer 11 flyby of Jupiter in 
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December 1974. The non-uniqueness of derived models will be examined 
utilizing the generalized inverse techniques that have been developed during 
the last decade in other branches of geophysics, notably, seismic data 
inversion (Wiggins, 1972), gravity and magnetics (Pedersen, 1977)* and studies 
of the electrical conductivity of the earth (Johansen, 1977). The generalized 
inverse techniques are of wide applicability in studying the often peculiar 
properties of linear (or linearized) systems and provide a framework for the 
evaluation of various models. The method underscores the uniqueness problems 
associated with the derivation of spherical harmonic expansions for the 
internal field and facilitates a systematic study of model parameter 
resolution. As the inevitable comparisons between models derived from 
different spacecraft flybys are made (e.g.. Hide and Malin* 1979), it is 
important to evaluate the observations for consistency and develop models 
which are applicable to all of the flybys. The ultimate goal is to find one 
model, or field description, which is appropriate to all the observations, 
allowing integration cf and sensible comparisons among the individual flybys. 
Thus, we will in addition examine the basic models chosen to represent the 
data and the validity of assunptions required by the conventional spherical 
harmonic analysis of these data. In this regard we will assess the importance 
of the magnetic field contribution due to external currents in the Jovian 
magnetodisc and propose improved models which exolioitly incorporate these 
effects. For this purpose, models of the inner magnetosphere resulting from 
an analysis of Voyager 1, Voyager 2, and Pioneer 10 observations (Connerney et 
al., 1981) will be used to supplement the conventional spherical harmonic 
analysis . 

The primary purpose of an internal field model is to provide an accurate 
description of the magnetic field throughout all space, not just along a 
single trajectory. Thus the ultimate test of a derived field model is how 
well it predicts the field in regions of space far removed from the locus of 
observations to which the model was fitted. Correlative data can be used to 
test various models, as well as additional in situ magnetic field 
observations. Comparison of predicted charge particle satellite absorption 
with observations (Acuna and Ness, 1976) is a powerful test of field models, 
since the trapped charged particle population reflects a global magnetic field 
topology rather than a localized observation. Indeed, satellite charged 
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particle absorption signatures can be used (Acuna et al., 1980) to constrain 
the initial field models. Other field topology diagnostics which have been 
applied to the Jovian field with varying success include ultraviolet auroral 
observations (Broadfoot et al., 1980), the polarization of Jovian decimetric 
radiation (Birmingham, 1980; de Pater, 1981) and the frequency and beaming 
pattern of decameter radio emission (Alexander et al., 1975). These 
observations can also in principle be used to constrain ix>ssible field models. 

METHOD OF ANALYSIS 

Increasingly more sophisticated models of the int<'-nal Jovian magnetic 
field have been proposed as the quantity and quality of available data 
increased. The Pioneer 10 data were initially modeled fis a centered, tilted 
dipole and an offset, tilted dipole (Smith et al., 1974; Van Allen et al., 
1974). Ultimately, data were interpreted utilizing the spherical harmonic 
analysis (Chapman and Bartels, 1940) conventionally used in studies of the 
earth's magnetic field. In the usual spherical harmonic analysis it is 
assuned that the data were obtained in a source free region of space (such 
that s 2 X B s 0). The magnetic field B can then be expressed as the 
gradient of a scalar potential function 


B s -7V = -7(V® + V^) (1) 

which can be written as a sum of two potentials, representing sources internal 
and external to the region of interest. The potential is expanded in terms of 
spherical harmonics 

V = V® + = a ! (r/a)" T^® + (a/r)"'"^ T^^ (2) 

n=1 


ivhere r is the distance to the planet's center, a is the planetary radius, and 
e i 

the T ^ and T^ are given by 

Tn^ = 2 p^'"(cos6) cos mij) + h^^'^ sin m<j»] (3) 

m=0 
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( 4 ) 


T * s 2 P_'^(cos6) CG„'" cos H ™ sin mi^] 

n n n ^ n 

msO 

where 9 and »? ars co-latitude and east longitude, respectively, the are 
the associated Legendre functions with Schmidt normalization, and the 
hn*”, are the internal and external Schmidt coefficients. The series 

(2) is terminated at a suitable n = and the coefficients h^"*, G^'", 

are chcsen to minimize, in a least squares sense, the differences between 
the model field and the observations. This is accomplished by solution of an 
overdetermined linear system 

Y = A X (5) 

where ? is a column matrix of the N magnetic field observations, S is an N x M 
matrix relating the observations to the model parameters, X, arranged as a 
length M colunn matrix of the g^”", h^™, G^”*, coefficients. The details of 
the usual least squares procedure in this application are summarized in Acuna 
and Ness (1976b). 

The shortcomings of this usual spherical harmonic analysis as applied to 
the magnetic field observations at Jupiter have been discussed by both teams 
of investigators involved (Davis and Smith, 1976; Acuna and Ness, 1976b). 
Briefly, these shortcomings fall into two distinct areas which will be 
discussed separately. The first involves inadequacies of the basic model and 
the second involves the often peculiar properties of linear systems. 

EXTERNAL FIELDS 

Inadequacies in the basic model include the assumption that the data were 
obtained in a source- free space and the inefficiency of modeling fields due to 
external currents as an expansion of the form (2) with a limited number of 
terms. In the spherical harmonic analysis of Voyager 1 and 2 magnetic field 
data (Ness et al . , 1979a, b) the reduced g^° terras were attributed to the 
repeated passage of these spacecraft through an annular current sheet (i.e., 
the magnetodisc) circling Jupiter. It is now clear (Connerney et al., 1981) 
that Pioneer 10 and Voyagers 1 and 2 were repeatedly immersed in a current- 
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carrying region even as they traversed the innermost Jovian magnetosphere. 

6 

Azimuthal current densi'-ies in the magnetic equatorial plane of 5 x 10 
2 

A/Rj at a radial distance of 5 Rj were typical of models considered for both 
Voyager and Pioneer 10 encounters. Consequently, it is necessary to analyze 
these data in the context of a model which is applicable in a region of space 
containing sources. While it is quite possible that Pioneer 11 remained in 
source-free apace during encounter, it is most probable that the magnetic 
field observations include a sizeable contribution from the external current 
system. Some of the ’systematic effects' in the residual (observed minus 
calculated field) noted by Davis and Smith (1976) in their analysis of the 
Pioneer 11 vector helium magnetometer data attest to the presence of this 
sheet perturbation field. For this encounter ti.e problem is one of 'leakage' 
or aliasing of the field due to external currents into coefficients describing 
the internal field (g„”i h^'”) . 

A reasonable starting point in the effort to construct Improved models is 
to concatenate a spherical harmonic expansion of the type (2) containing only 
internal terms with an explicit model of the current system in the Jovian 
magnetosphere. While several models of the Jovian magnetosphere and external 
current system exist (Barish and Smith, 1975; Goertz et al., 1976; Jones et 
al . (1980); Beard and Jackson, 1976; Engle and Beard, 1980; Connerney et al., 
1981) only the model of Connerney et al. (1981) represents a detailed vector 
fit to magnetometer observations in the innermost magnetosphere (< 20 Rj) , the 
region of importance here. The other models represent qualitative fits (Beard 
and Jackson, 1976; Engle and Beard, 1980; Barish and Smith, 1976) or are 
applicable in the more distant magnetosphere (Goertz et al., 1976) and often 
match only scalar magnitude. The Euler potential models of Goertz et al . 
(1976) and Jones et al . (1980) are valid only in the near equatorial region of 
the distant (R > 20) magnetosphere, while the Biot-Savart model of Jones et 
al . ( 1980) was compared only with observations obtained at radial distances 
exceeding 20 Rj. Thus the Connerney model , particularly the analytical 
approximations discussed in the Appendix to Connerney et al. (1981), will be 
utilized to describe the external field. In a combined model we write 

B s B^ + b 
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where B is the total field, B the perturbation field due to external currents, 
and ia the internal field. ia derivable from a scalar potential (1) and 
B ia computed as the curl of a vector potential 5 due entirely to external 
currents, ITie parameters of the model external current system and model 
internal field can then be simultaneously determined by inversion of the 
magnetic field data. In this paper we treat the parameters of the modeled 
external current system as fixed constants, determined by the VI, V2, and P10 
encounters. Thus we assune at present that the near axis field due to 
external currents in the Jovian magnetosphere is time invariant, such that the 
model field fitted to the VI, V2 and P10 observations is appropriate to the 
P11 flyby as well. The magnetic field of the magnetopause and tail current 
systems is estimated to be only a few gammas in magnitude in ttie inner 
magnetosphere (Engle and Beard, 1980; Ness et al., 1979c) and is neglected in 
the computation of b. 


The perturbation field, b, is computed using the approximate formulas 
given in the appendix to Connerney et al , (1981) for the near-axis region of 
the Jovian magnetosphere (in magnetic equatorial coordinates): 


>* 0^0 

B = ( ) n/F. - l/F,] 

^ 2 2 ' 


and 


B, = ( °?) C2D (Z‘ + a‘)‘ 


4 F ^ f' ^ 

^ M *^2 


where 


F^ = [(z-D)^ + 


Fg = C(z+D)^ + 

p and z are the radial and vertical coordinates, D is the annular current 
sheet half-thickness , and a (b) is the inner (outer) radial extent of the 
current annulus. For the VI sheet model used in this paper, we adopt 
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1 225, as5Rj, bsSORj, and D = 2.5 Rj, the resulting field quantities 
expressed in gammas. Magnetic field quantities in the usual System III (1965) 
coordinate system are obtained by the appropriate coordinate transformation. 
The magnetic field topology in the magnetic equatorial plane of the field due 
to the VI model current is illustrated in Figure 1. - The magnitude of the 
field due to the sheet currents at the origin is -r 200 y. 

GENERALIZED INVERSE THEORY 

The analysis of a linear system of the type (5) is conveniently 
accomplished using the generalized inverse theory that has been successfully 
utilized in the last decade in a nimber cf branches of geophysics (Jackson, 
1972; Pedersen, 1977; Wiggins, 1972). Instead of proceeding directly with the 
linear system (5), however, we will present the methodology in a form 
applicable to non-linear problems, in anticipation of non-linear relationships 
between the observations and the model parameters. For example, field 
magnitude observations (obtained while the spacecraft traversed Jupiter's 
shadow and lost orientation information) are non-linear in the parameters of 
the usual spherical harmonic models. Vector field observations are 
non-linearly related to model parameters in the context of the kind of models 
suggested in the previous section. Indeed, should correlative observations be 
included in the inversion process, it is likely that they too will be 
non-linear in the model parameters. 

We assume that the i-th observation, Y^_ , is related to the model 
parameters by the function F’j^(Xj), The functions ^^^(Xj) may be expanded in a 
Taylor series about some initial parameter set, Xj°, 

0 

Y, = F,(X °) + _| 4X, + ... 

J 

Neglecting all but terms of order 1, and letting 
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( 5 ) 






3X, 


AX, 




which has the sariie form as the matrix equation 


7 . S 


If we let the oolunn vector 7 represent the model residuals (observed minus 
modeled field) , the vector X represent the parameter corrections required to 
bring the model into closer agreement with the data, and the 5 matrix, a 
matrix of partial derivatives of the model field with respect to the model 
parameters, i.e., 








3X^ 



We will proceed with (5) with the understanding that non-linear problems are 
accommodated by considering them to be locally linear and iterating to a final 
solution . 


The matrix formulation of the generalized inverse method utilizes the 
singular value decomposition of Lanezos (LanezoSy 1961) to rewrite (5) as 

7 = U A X (6) 

where 0 ig an N x M matrix consisting of. the M orthonormalized eigenvectors 

mmmmT 

associated with the M largest eigenvalues of AA , 7 is the H x M matrix 

T — 

consisting of the orthonormalized eigenvectors of 3 as columns, and a is an 

—T •• 

M X M diagonal matrix consisting of the eigenvalues, of A 3. The a matrix 
is by convention assembled with the largest X^ in the upper left, all elements 
positive nqd in order X,j > Ag > ... Xj^, i.e.. 
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A = O Xj 
T 

Opdrtting on the left with d , we write 

(7) 

Since A is an M X M diagonal matrix (7) can be regarded as M independent 
equations relating 'eigendata' (on the left), through the eigenvalues Xj ♦ to 
eigenvectors of parameter space, the linear combination of the original 
parameters The solution to (5)» that is, the parameter vector X 

minimizing in a least squares sense (Lanozos, 1961) the differences between 
the model and observations, is given by 


= 7 A""' ? (8) 

* T 

writing 3 s d 7, the solution can be constructed by a summation over the 
orthonormalized of parameter space: 

M 

X = Z (S^/\^) Vj_ (9) 

i*1 


Considerable insight into the estimation problem can be obtained by an 
examination of (9). Assune that in the original statement of the problem 
(equation 5) the data are (by suitable transformation — see Jackson, 1972) 
normal random variables of zero mean and unit variance. In (9) the are 
then of unit variance. The variance of the parameter vector (solution) is 
largely due to just a few of the x^*s of the A matrix corresponding to the 
directions in the parameter space (V^’s) along which the solution is poorly 
determined. The singular value decomposition yields not only a method of 
constructing solutions but a characterization of both parameter and data space 
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as well. 


V/e stated earlier that magnetic field observations along a flyby trajec- 
tory are insensitive to certain linear combindtion of parameters. One 
advantage of the singular value decomposition is that these parameter vectors 
which are poorly constrained by the available observations are explicitly 
identified, i.e., they are the eigenvectors associated with the small 
eigenvalues of ff K, jo illustrate this, we show in Figure 2a a surface 
isointensity contour map of a maliciously constructed, hypothetical Jovian 
magnetic field. Illustrated in Figure 2a is the surface expression of a 
particular combination of internal field parameters (eigenvector of parameter 
space) which leads to a virtually unobservable magnetic field along the 
Pioneer 11 trajectory, but a significant field elsewhere. Had Jupiter's 
internal field been so constructed, the Pioneer 11 fluxgate magnetometer would 
not have detected it. We have in this example assutsed observations within 5 
Rj at .5 minute intervals and a random noise component of .005 G everywhere; 
the parameter set for this example is just the appropriately scaled 
eigenvector associated with the minimun eigenvalue. In Figure 2b we show the 
surface isointensity contour map of an internal field which would be 
undetectable assuming measurements with a random noise component of 15 of the 
ambient Jovian field magnitude rather than a constant noise component as in 
2a. Observations within 5 Rj of Jupiter at ''' 1 minute intervals are assumed. 
Note that in this example the field magnitude is smaller but comparable to the 
previous example, and that the eigenvector is considerably different , 
reflecting a difference in information content in the two data sets. 

The A matrix is a function of the trajectory, the model, and the 
transformation utilized to condition the data. If the model chosen is 
sufficiently accurate, or descriptive of the magnetic field environment, the 
(observed minus modeled field) will be of zero mean. Each observation, , 
and correspondingly, each row of the S matrix is divided by the estimated 
standard deviation of that observation, such that the standardized 
variable Y^^/o^ has a unit standard deviation. Thus the residual for each 
observation is scaled with the width of the population from which it was 
drawn, its expected error. For independent observations, the A matrix assunes 
the form 
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( 10 ) 


1 3F- 1 3F, 

A' = ' 



1 — 




corresponding to a least squares minimization of 

(Y - A 5 (Y - A 


( 11 ) 


v/here the matrix C is the diagonal matrix with elanents s 1/<J^^ 


For data vAiich are not statistically independent, the inverse covariance 
matrix of the observations replaces 5 in (11) (Jackson, 1972). In either case 
the choice of weights relates to the statistics of the noise on the data, and, 
by modification of the A matrix (10), alters not only the final model but the 
orthogonal Ization of data and parameter space ast well. Thus care must be 
taken to insure chat the weights chosen reflect some real knowledge of the 
noise statistics of the data. 

_T- 

We have thus far assumed t.hat the M eigenvalues, of the matrix A A 
are non-zero, in which case (9) represents the classical least squares 
solution. The Lanczos inverse also leads to a solution when only K of the M 
eigenvalues are non-zero. The non-zero are arranged in the diagonal matrix 
A as before with x^^, i > K = 0. Equation (9) becomes 

K 

X (K) = Z (12) 

i=l 


i.e., only the eigenvectors corresponding to non-zero eigenvalues are allowed 
in the solution. Those parameter combinations which, by (6), do not 
contribute to the observations are ignored in the construction of the solution. 
The minimum length solution (12) is unique in the space spanned by the K 
eigenvectors , i s 1 , K, and has no projection along the remaining M-K 
eigenvectors (corresponding to the zero eigenvalues) which define the region 
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of parameter space that is inaccessible to the present linear system. This 
inaccessibility reflectes a basic lack of information inhererst in the 
available observations. The unavoidable consequence of such a deficiency is a 
lack of parameter resolution, that is, the existence of linear combinations of 
parameters (Vj^, is K+1 ... M) which may be added to the solution of the 
linear system without change to the model response. 

Eigenvalues which are non-zero but small lead to a similar loss of 
resolution when noise on the data is considered. The projection of the 
solution vector onto the corresponding to the large Xj, will be well 
constrained; the projection onto the associated with the small will be 
relatively poorly constrained, by virtue of the factor (fl^/X^) in (12). 

Errors in the parameters resul'.ting from noise upon the observations will 
appear predoroinently in linear combinations corresponding to the poorly 
determined eigenvectors. A vivid and perhaps more familiar example of this 
effect occurs in the geomagnetic field models derived from scalar 
observations. Stern and Bredekamp (1975) have shown that errors in field 
models derived from scalar data are enhanced for certain sequences of terms, 
first derived by Backus (1970). This 'Backus effect', subsequently observed 
(Stern et al., 1981) in models derived from Magsat observations, suggests the 
near- singular nature of the A matrix for this problem, and the presence of 
eigenvectors (related to the Backus series) which have small associated 
eigenvalues. Again assuning the data have been transformed such that the 
are statistically ir, dependent and of unit variance, the variance of the jth 
model parameter is 


K<M 

S^. (K) >. Z 



( 13 ) 


In constructing a solution vector, we are free to restrict the solution to K < 
M eigenvectors, reducing the parameter variance to an acceptable level. In so 
doing one avoids a large contribution to the parameters due entirely to a 
small noise component on the data, at the expense of an unavoidable loss of 
model parameter resolution. 
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Paraneter resolution is conveniently described by the resolution matrix, 
R, (e.g., Jackson, 1972) relating the K eigenvector parameter estimate )J(K) to 
a solution of the linear system (5), 

!») = 7 An"' ij'' r 

s 7 Ajj"’ 7^ 0 A 7'^ ii 
. 7, 7j ^ 

The subscript K to a matrix denotes the matrix obtained by setting each colunn 
i for i > K of the original matrix to zero. An element Xj^(K) of the estimated 
solution is the convolution of the ith row of the resolution matrix with a 
solution of (5). For K s M, the R matrix is the identity matrix; as fewer 
eigenvectors are admitted in the construction of the solution, the R matrix 
off-diagonal elements grow at the expense of the diagonal elements, reflecting 
a loss of parameter resolution. An example illustrating the trade-off between 
parameter standard deviation and resolution is given in Figure 3. In this 
example we assume an internal spherical harmonic field model of order 3 and 
500 vector magnetic field observations along the Pioneer 11 trajectory within 
5 Rj of Jupiter, the observations including a noise component with .005 G 
standard deviation. The RMS model residuals... g.|^ parameter resolution 
(diagonal element of R) and the g.^° parameter standard deviation are shown as 
a function of the number of eigenvectors admitted in the construction of a 
solution. The very modest improvement in model fit (RMS) attained by using 
the last two eigenvectors is accompanied by a large increase in parameter 
standard deviation and only a small gain (25) in parameter resolution. An 
estimate for gy of 4.35 G, based on the 13 ev fit, may be more appropriate 
than the estimate of 4.18 G resulting from the 15 ev solution. With 13 ev , 
the off-diagonal terms of the resolution matrix are sufficiently small that 
'leakage* frcm the higher order parameters (which are presumably small in 
magnitude) inco the g.^^ estimate should not be important, i.e., is 
expected to be adequately resolved at R^.j = .98. 
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The advantages of obtaining additional magnetic field information far 
removed, from the observations available along the flyby trajectory can be 
understood by considering only the poorly determined parameter vectors. Such 
information may include the correlative observations referred to earlier as 
well as additional in situ magnetic field observations. The additional 
observations will in general be in regions where the poorly determined 
eignevectors lead to relatively large fields (e.g., Figure 2). These new 
observations effectively limit the range accessible to the previously poorly 
determined eigenvectors, acting as a constraint on the solutions. This can be 
quantitatively explored by forming a combined system of all the observations 
(e.g., equation 10) and comparing the new eignevalues with the old, or in an 
approximate way (Jackson, 1972) by estimating the (larger) eigenvalues of an 
augmented A raatriy, of N 1 rows. The potential benefits in terms of improved 
model parameter variance and resolution are great since they depend critically 
upon the minimum eigenvalue. 

To define the entire range of values a parameter may assume, we use the 
concept of an extreme parameter set as developed by Johansen (1977). Because 
of correlations between the parameters, the extreme value of any parameter is 
achieved for a specific parameter vector — its extreme parameter set. This 
represents he direction in parameter space along which the model response 
(RMS change) is minimized for a change m a specified parameter. This 
direction is found by considering the change in the model response, ^Q, due 
to an excurcion, e, in parameter space from the model minimum 

aQ = V f- e (14) 

In the coordinate system defined by the eigenvectors 7^, AQ is given by 

AQ = ? (15) 

describing an M dimensional hyperellipsoid with axes parallel to the 7^'s of 
length a 1/X^. The extreme parameter set associated with the extreme values 
of the ith parameter is determined by requiring the gradient of the error 
surface to have only a component along the ith parameter. 
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7 = sa^ 

where is the column vector with the ith element equal to one and all others 
zero, resulting in 


Ej^ s 3 V a”^ 7’^ (16) 

where S is a constant and E^^ is the extreme parameter set for the ith 
parameter. Any linear combination of the model parameters may be extremized 
in this fashion. Zn actual practice one perturbs the solution in the 
direction given by (16) until Ag appropriate to a specified significance level 
is achieved to find the magnitude of the extreme parameter sets, as suggested 
by Johansen (1977). This procedure requires that the problem be sufficiently 
linear in the neighborhood of the solution such that the basis vectors of 
parameter space (V^'s) are not appreciably altered in the process. The 
problems considered hcfre are sufficently linear that we may use for S the 
value obtained from the linear approximation, resulting in 


E^ s V a“^ f a^/ll a"'' f II (17) 

for AQ s 1, corresponding to a 685 confidence limit. 

One additional practical consideration relating to the non-linearity of 
the problem has to do with the method (9) of constructing solutions at each 
iteration. If the A matrix is close to being ill-conditioned in the sense 
that some of the are very small, the solution vector will require a large 
step in parameter space which may well be greater than the region in which the 
linearization is appropriate. The iterative process may then diverge unless 
some method of limiting the iterative step size is employed. In such cases 
the stability of the iterative procedure is improved by limiting the step size 
(noteably in the directions corresponding to the small eigenvalues) by using 
(e.g., Lawson and Hanson, 1974, Chapter 25) 

K ^i H - 

X = r V (18) 

i-l + a^) 
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where a is chosen approximately equal to the snallest eigenvalue included in 
the inversion. This technique (known as Marquardt’s algorithm) is invoked 
only as convergence problems arise and is quite unnecessary in the context of 
the present nearly-1 inear problems considered here. 

APPLICATION TO PIONEER 11 OBSERVATIONS 

The Pioneer 11 spacecraft had onboard two magnetic field experiments 
(Smith et al., 1976; Acuna and Ness, 1976), each of which recorded the 
magnetic field along the Pioneer 11 trajectory. The vector helium 
magnetometer (VHM) has been described by Smith et al. (1975) as accurate at 
the 1% level (at low field values) whereas the high field fluxgate (FGM) is 
essentially a constant noise instrument with a quantization uncertainty of 
.006 gauss (Acuna and Ness, 1975). The high field inctrunent would be 
expected to provide a better estimate of the Jovian field near close approach 
(jBl V' 1.13 Gauss at the inbound ocoulation) while the vector helium 
magnetometer would out perform the fluxgate in a low field environment. In 
either case the random noise component assumed in what follows is a sum of all 
noise sources, including, for example, instrumentation noise and local 
magneto spheric noise. A single least squares fit to the combined data sets, 
each observation appropriately weighted as discussed in the previous section, 
would result 'in the optimal model. Unfortunately, the t\x> data sets are 
incompatible in the sense that the systematic differences in the observed 
field as measured by the two instrunents are larger than expected on the basis 
of the instnment descriptions. This is illustrated in Figure 4, which shows 
the differences between the VHM and FGM measurement, for each component of the 
field during encounter as a percentage of the total field magnitude. Rather 
than produce a model which fits neither set of data, we will of necessity 
continue the practice of considering each data set independently. Resolution 
of the differences between the two data sets is beyond the scope of the 
present work. 

The FGM observations have been analyzed by Acuna and Ness (1976a,b) in 
the traditional least squares fit to a spherical harmonic expansion of order 3 
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including internal terras only. They obtained the GSFC 0^ internal field model 
as the unweighted least squares solution to 685 vector observations at radial 
distances of R < 6 Rj. Details of their analysis are given in Acuna and Ness 
(1976b). We will first present the results of a similar analysis, assuning an 
expansion of internal harmonics to order 3, but within the framework of the 
generalized inverse methodology. A model explicitly including the effects of 
the field due to external currents will then be considered. 

Our FGM data set consists of 499 vector observations obtained 
approximately every 35 seconds along the Pioneer 11 trajectory within 5 Rj 
radial distance, and 20 field magnitude observations obtained near closest 
approach vihen the orientation of the spacecraft is uncertain. Inclusion of 
the magnitude observations along this part of the trajectory increases the 
Information content of the data set at the expense of introducing some 
non-linearity in the estimation problem and requiring an iterative solution. 
Since the FGM is essentially a constant noise instrunent (Acuna and Ness, 
1975), and the FGM observations are limited by instrunent noise, we choose for 
each observation equal weights (Eqn. 10) of s ,005 Gauss. The use of equal 
weights introduces only a scale factor as compared with an unweighted least 
squares and does riot effect the orthogonallzation of parameter space or the 
final model. Assuning as a model a Schmidt normalized spherical harmonic 
expansion of order n = 3 (internal terras only) we find after two iterations a 
15 eigenvector fit (Eq. 9) very similar to the GSFC 0^^ model, differences 
between the two solutions due only to a difference in observations used. The 
salient features of the singular value decomposition are illustrated in Figure 
5, at the final iteration. In this figure the eigenvalues of the A matrix, 
and the associated eigenvectors of parameter space are listed, in addition to 
the solution vector, parameter standard deviations (Eq. 13) and extreme 
parameter sets (Eq. 16). A diagnostic parameter 'lin' is also indicated: this 
is the cosine of the angle between the parameter vector and the local normal 
to the error surface. A change in the value of an individual parameter with 
lin “ 1 results in a comparatively large model response (RMS increase) since 
that parameter vector is nearly perpendicular to the error surface, that is, 
it is oriented in the direction of the maximum model response (RMS change). 

The weighted RMS residual is computed at each iteration assuning internal 
consistency (Birge, 1932) 
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essentially a prediction of the probable error based on propagation of errors, 
and also assuming external consistency 
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based on the realized differences, c^. A comparison of these two estimates is 
often useful in identifying the presence of systematic effects in the data. 


The eigenvalues for this system range from <^50 to The classical 

measure of the condition, or stability, of a linear system is the 'condition 
nunber' , v, defined as the ratio of the largest and smallest eigenvalue 
( Lane zos , 1971), 


or sometimes the square root of the above. In the present example, a 
'condition nunber' of v * 65 is interpreted in the following way: errors in 
the 15th generalized parameter (15th eigenvector) can be expected to be v/'65 
times larger in magnitude than errors in the 1st generalized parameter. The 
'condition nunber' is thus a useful tool in diagnosing the need for a singular 
value analysis, but it cannot address individual parameter errors or 
resolution, or correlations among the model parameters. 

Inspection of the eignevectors associated with the eigenvalues of the A 
matrix is useful in visualizing the characteristics of the problem solution. 
For example, in Figure 5, the eigenvectors associate with the small eigen- 
values are the linear parameter combinations which are poorly constrained by 
the observations. The parameter with the smallest absolute error (J =2; 
parameter g^ ) has only small components among the last (poorly determined) 5 
eigenvectors, that is, it is well approximated by a linear combination of the 
first (well determined) 10 eigenvectors. Conversely, the parameters with the 
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2 1 

largest absolute errors (J = 1 1 , 82 . and J = 13i h 2 ) are those which are 
heavily represented in the last few eigenvectors. Furthermore, these two 
parameters are equally represented in the 15th eigenvector but with opposite 

signs, indicating a strong negative correlation between these two parameters. 

2 

A change in the parameter g- will most likely be accompanied by a nearly 

equal and opposite change in the parameter • A similar relationship exists 

1 3 

between the parameters g^ (J * 10) and g^ (j s 12 ), as evidenced in the I4th 
eigenvector . 

The extreme parameter sets are listed in Figure 5 for each parameter, 
computed using (16). These are the parameter combinations which may be added 
to or subtracted from the optimal solution to achieve the extreme values of 
each individual parameter. The error surface chosen here is AQ = 1 
corresponding to a 68% confidence limit. 

A useful illustration for the purpose of emphasizing the differences 
between internal field models is a surface isointensity contour map such as 
that in Figure 2a. Figure 2a, the ’undetectable planet’ was constructed by 
selecting the internal field coefficients to be proportional to the 15th 
eigenvector, that combination of parameters to which the observations are 
least sensitive. It is perhaps more instructive to compare the surface 
isointensity contour map corresponding to the best fitting model (Figure 6) 
with the surface isointensity contour maps that result when this parameter 
vector is added to (Figure 7a) or subtracted from (Figure 7b) the least 
•squares solution. These two extreme models lead to differences in the surface 
field magnitude of >/'1 G. 

While this example is representative of the differences in field 
topology and magnitude appropriate to the range of models consistent with the 
observations, it is by no means unique. Any linear combination of the poorly 
determined eigenvectors, consistent with the constraint (14) and some chosen 
significance level (AQ), can be added to the least squares solution with 
similar results. The surface isointensity contour map corresponding to the 
least squares 15 eigenvector solution,, is practically identical to that of the 
GSFC model of Acuna and Ness (1976), reflecting only differences in the 
radial extent of observations included in the analysis and the inclusion of 



additional magnitude obssrvations near close approach. Thus, the model 
uncertainties illustrated in Figure 7 can also be regarded as typical of the 
0^ model. 

For some applications a model constructed of only 13 eigenvectors, 
ignoring the last two (poorly determined) eigenvectors, may be preferable to 
the 15 eigenvector solution. In Figure 3 it is clear that the 15 eigenvector 
fit is only marginally better than the 13 eigenvector fit in terms of the RMS 
of the residuals (although the I4th and 15th parameter vectors are statisti- 
cally significant, in that 0^^ and 0^^ > i). The loss in parameter resolution 
that results in constructing the solution with only the first 13 eigenvectors 
is illustrated in Figure 8, which lists the resolution matrix as well as the 
quantities discussed in the context of the 15 eigenvector fit. All of the low 
order terms (ns 1, n s 2) are well resolved, since the corresponding diagonal 
elements of the resolution matrix are near unity. Additionally, the g^° 

(J s 9) is well resolved, and the term (J = 15) somewhat leas so. The 
remaining 3rd order terms are rather poorly resolved. The surface 
isointensity map of this model field, illustrated in Figure 9, is more dipolar 
than that of the 15 eigenvector fit as a result of the concentration of higher 
order paraneters in the unused eigenvectors. While the 13 eigenvector model 
may provide a better estimate of the well resolved parameters, it is unique 
only in a subset of parameter space. We are free to add linear combinations 
of the excluded eigenvectors to the solution, subject to some maximun 
allowable RMS increase corresponding to a chosen significance level. 

Two observations remaining to the model residuals for the FGM data set 
are of importance in regard to the significance level corresponding to the 
choice of AQ s 1. The first is that the RMS of the residuals, >^500 Y, is 
largely due to the quantization step size of the instrument; the true 
instrumentation noise is a fraction of that. The second is that the 
quantization noise is not random, rather, it exhibits a positive correlation 
at small lags and a negative correlation at larger lags. The correlation 
length is several samples (minutes) at large radial distances (5 Rj) where the 
field is comparatively small. At small radial distances the correlation 
vanishes. The quantization errors, however, do not seriously affect the 
model, as they are as often positive as negative. In constructing solutions 
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(9) only the quantity a * 0 7 is important, and the sum over the quantization 
errors is likely to be small since each positive residual is paired with a 
nearly equal negative residual in the immediate vicinity. The quantization 
errors contribute substantially to the RMS of the residuals, but not 
appreciably to the model misfit. Thus the correlation among the residuals, 
effectively dwreasing the nunber of independent observations, is largely 
offset by the decrease in the estimated noise upon the data (disregarding the 
quantization 'noise'). For this reason we regard the choice of o = 500 y, 
corresponding to 6855 confidenee level of a Gaussian distribution as a 
conservative estimate. 

We have demonstrated that a sizeable parameter vector can be found (e.g.. 
Figs. 2a, 2b) which leads to a small field along the Pioneer 11 trajectory. 

In much the same way a small field (systematic 'error') along the Pioneer 11 
trajectory can lead to a large error in the estimated parameters. The 
magnetic field of the large scale azimuthal currents produce a large ('/'200y) 
correlated 'error' along the entire Pioneer 11 trajectory under consideration 
(Connc.-ney et al., 1981). This 'error' contributes substantially to 8 = 0^7 
due to the large positive correlation (no cancellation) and thus results in a 
significant bias in the derived internal field models. The bes^ way to remove 
such a bias is to include explicitly the magnetic field produced by these 
currents in the physical model chosen to represent the observations. A 15 
eigenvector fit to the FGM observations, assuning a N = 3 internal spherical 
harmonic expansion and the VI model external field, results in a modest RMS 
decrease to 492 y but a significant change in many of the model paraneters 
(Table 1). Comparison of the surface isointensity contour map of this model 
field (Figure 10) with that of the earlier 15 eigenvector model (Figure 6) 
indicates surface field magnitude errors of «/'1 Gauss result from neglect of 
the external field in the earlier model. Differences in the internal field 
parameters resulting from inclusion of the sheet field are as large as .12 G 
(Table 1) and often larger than the corresponding parameter standard 
deviations. Thus the presence of an unmodeled field of external origin is at 
least as important as the noise on the data in this examaple. Since no 
additional free parameters have been introduced (the parameters of the current 
sheet regarded as fixed, determined by the more advantageous Pioneer 10 and 
Voyager flybys) the singular value decomposition remains unchanged. A 
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formulition of this problem In which the sheet parameters are treated as 
variables would show more clearly the interdependence of the internal 
coefficients and the current sheet parameters. 

The VHM observations have been analyzed (Smith et al., 1976; Davis and 
Smith. 1976; see also Smith and Gulkis. 1979) in the traditional least squares 
fit to a spherical harmonic expansion. Noteable differences in their 
analysis induce a larger radial range of observations used (r < 8 Rj), the 
use of external spherical harmonica in the expansion (2), and the use of a r* 
weighting scheme whereby all observations are weighted with some power of the 
radial distance of the observation. For reasons discussed in Davis and Smith 
(1976) this group prefers a weighting of each observation. In addition, 
they use (Davis and Smith, 1976) minute averaged values of the observed 
field along the trajectory in their analysis. 

Our VHM data set consists of 324 vector (</’ 1 minute averages) 

observations within R « 5 Rj of Jupiter as obtained from the National Space 

Science Data Center at Goddard Space Flight Center. In contrast to the 

uniform weighting appropriate to the FGM observations, we adopt for the VHM 

observations weights proportional to the local field magnitude, = .01 jBj. 

The constant of proportionality adopted here (.01) enters only as a scale 

factor in the error analysis. This particular weighting, appropriate for 

observations with a noise component of I* of the local iB| , results in a more 

heavily weighted (in r) system than any of the r* schemes proposed by Davis 

and Smith (1976). (The weights, w^, defined by Smith et al . (1976) and Davis 

and Smith (1976) appear linearly in their formulation of the least squares 

problem. Therefore these w^ must be compared to the square of the weights 

(1/a^) appearing in this work.) The resulting model (listed in Table 1) leads 

to a surface isointensity contour map, illustrated in Figure 11, that is much 

more similar to that of the FGM (or GSFC Om) model than previous analyses of 

3/2 

these data (using r weights) would suggest. The choice of weights 
proportional to ! B) affects not only the final model but the orthogonal ization 
of parameter space as is illustrated in Figure 12. We have assumed in this 
example a model field of internal terms to order 3 so that the VHM and FGM 
inversions can be directly compared. Comparing the eigenvalues of the two 
inversions, we note that the FGM observations lead to a slightly better 
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determination of some eigenveotoriJ and a poorer determination of others. Yet 
individual parameters would be better determined by the VHM inversion since 
they require a summation over many eigenvectors. The 'lin' diagnostic 
parameter demonstrates that all of the dipole parameters i;i the VHM inversion 
lie closer to the eigenvectors of parameter space in this problem and 
therefore are more easily estimated independently of the other parameters, ^ 
detailed study of the remaining features of this singular value decomposition 
is similar to the PGM example and thus will not be repeated; instead we will 
address an additional concern of systematic errors present in the VHM model 
residuals. 

The fact that the model residuals are characterized by an extremely large 
correlation length precludes any attempt to quantify the parameter errors. 

For this reason the extreme parameter sets have been omitted and the quoted 
parameter standard deviations are to be disregarded. These systeriatic errors 
include the presence of an unmodeled external field, which we shall return to 
shortly, and systematic errors associated with instrumental range changing and 
other effects (see Davis and Smith, 1976, Figure 1 and accompaning 
discussion). While these effects are only as large as 'J' .5 or IS of the total 
field they do dominate the observed residuals. Another systematic effect in 
previous analyses of these data involves the use of 5 minute averages of the 
field as instanteous values. Near close approach, the 5 minute averaged field 
differs by as much as 0.6X from the instantaneous field at the midpoint of the 
data averaging interval (which is used as the model field). Since this effect 
is comparable in size to the observed residuals it seems appropriate to use 1 
minute averages instead or compute the model response as the averaged field 
over the measurement Interval. We have chosen to use >^1 minute averages as 
supplied by the NSSDC for computational convenience. 

Inclusion of the magnetic field contribution due to the external currents 
in the analysis of the VHM observations leads to a much improved fit to the 
observations. The resulting parameter set, listed in Table 1, leads to a 
surface isointensity contour map (Figure 13) that still resembles the previous 
map but differs in field magnitude by ^1 G at the surface. Differences in the 
internal field parameters are as large as ^ .2 G, occurring predominantly in 
the quadrapole and oetopole terms. A more dramatic decrease in the RMS of the 
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residuals from 93 y to 57 t results from the use of the VI current sheet model 
in the Inversion, in contrast to the relatively modest Improvement (50^< y to 
492 y) noted in the FGM fit. It is worthwhile to note again that this 
improvement in the RMS residuals is not gained at the expense of additional 
free parameters, i.e., the sheet field is regarded as a known (based on the 
other Pioneer 1& and Voyager 1, 2 flybys) quantity. A relatively modest 
further improvement in the RMS of the residuals (to 56y) is obtained by 
assuning either a closer inner radius of 4 Rj rather than 5 as in the VI 
model, or a constant » 250 rather than 225 as in the VI model. Again, 

a simultaneous Inversion of the internal and external field parameters would 
yield more information, but must await further analysis of the Voyager 
observations. The kind of correlations existing between internal and external 
field parameters that would be evident in a singular value decomposition of a 
combined internal/external system are already implied in the present analysis. 
Note that the internal field parameters have been able to ’absorb’ most of the 
external field, which averages >^200 y along the trajectory. The ability of 
the external field to partially masquerade as internal coefficients, leading 
to the large parameter changes indicated in Table 1, portend.*, such an 
interplay of internal and external field parameters. 

DISCUSSION 

The internal Jovian field models presented here differ substantially from 
the models of Acuna and Ness (1976) and Smith et al. (1976). In the first 
case the difference arises almost entirely out of a consideration of the 
fields of external origin, while in the -econd both the weights applied to the 
observations and the inclusion of the fields of external origin contribute to 
the difference. The resulting models are found to be more similar than 
previously thought. The large differences in the models of Table 1, based on 
the VHM observations, and the models of Smith et al. (1976) bears witness to 
the effects of the weighting scheme utilized to condition the observations , 
and to the increased importance of external fields in this data set. Indeed, 
the weights utilized in this pjiper as appropriate to the VHM observations have 
in fact made the model more sensitive to small fields at greater distances and 
elevated the importance of unmodeled external fields. 
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The differences between the internal field models advocated here (GSFC 15 
CVS end JPL 15 evs, outlined in Table 1) and previous models are more readily 
apparent in Figures 14 and 15. In Figure 14 we compare the GSFC 15 evs and 0^ 
models on the (oblate) Jovian surface by contouring the difference in field 
direction and the normalized difference in field magnitude resulting from each 
model. These differences are greatest in the northern hemisphere, where the 
local field angle diffen by up to 4° ^ind the local field magnitude differs by 
as much as In Figure 15 the same comparison is made between the JPf. 15 
evs model and the SHA 23 model. He>*e again the differences appear greatest in 
the northern hemisphere, where the local field direction differs by up to 30° 
and the local field magnitude differs by as much as 5011. A similar comparison 
between the JPL 15 evs and the more recent P11A model (Smith et al . , 1976) 
leads to differencs^ contours much like those illustrated in Figure 15 but with 
maximun differences of 18° in field angle and 20 % in field magnitude. The 
P11A model is a fit to observations obtained by both Pioneer 10 and 11, 
however, and therefore is not directly comparable to the JPL 15 evs model. 

But the inclusion of additional observations along another trajectory in the 
P11A model fit apparently results in an internal field model that more closely 
resembles the JPL 15 evs model presented here. 

The generalized inverse technique described here should prove valuable 
not only in the estimation of planetary raagneic fields and the study of 
non- uniqueness problems encountered in geomagnetic models based on scalar 
observations (e.g., Stern et al., 1981), but in a variety of other 
multivariate problems as well. The methodology has been widely used in 
related geophysical disciplines in the past decade and has proven to be a 
powerful analysis tool. Indeed, some (e.g., Wiggins, 1972) consider 
multivariate inversion studies without such resolution analyses to be 
incomplete at best, and often misleading. In such cases the magnitude of the 
condition number of the (appropriately weighted) matrix to be inverted can be 
diagnostic of the need for a singular value decomposition. 

The kind of model advocated here, combining an internal spherical 
harmonic expansion with an explicit model of the field due to external current 
systems, is regarded as essential to understanding and integrating the 
magnetic field observations of each of the spacecraft that have thus far 
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encountered Jupiter. Perhaps one model will ultimately enable meaningful 
comparisons of the various flybys and a more detailed description of the inner 
magnetosphere of Jupiter. Further progress can be expected from analysis of 
the Voyager observations and also possibly x'rom further analysis of the 
Pioneer 10 and 11 observations. At the present level of analysis, the small 
instrunental and other systematic effects present in the VHM observations are 
the most important remaining source of errors. For the purpose of eliminating 
or reducing these effects (and quantization errors in the FGM data as well) 
further analysis, conducted in payload coordinates, would be extremely useful. 
Analysis of the data in payload coordinates, while cunbersome, enables the 
modeling of quantization errors, axis non-orthogonality, and range sensitivity 
or gain errors. Such an analysis is not curr-ently possible with the reduced 
data available at the NSSDC. 
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FIGURE CAPTIONS 


Figure 1, Magnetic field topology in the magnetic equatorial plane of the 
field due to large scale azimuthal currents in the Jovian magnetosphere 
(VI model of Connerney et al., 1981). The current flows eastward in an 
annulus 5 Rj thick extending from 5 to 50 Rj. 

Figure 2. 'Invisible planet': Surface isointensity contour map of a Jovian 
model magnetic field which would not have been detected by a magnetometer 
on Pioneer 11 (a) (upper panel) assuning observations with a .005 G 
random noise component; (b) (lower panel) with a noise component 
proportional to the local field magnitude (see text). A dynamical 
flattening of 1/15.4 is assuned in the determination of the surface 
equipotential (all surface maps). 

Figure 3. Pioneer 11 FCM inversion example illustrating the trade-off between 
parameter resolution and parameter standard deviation of the main dipole 
coefficient (g^° term). 

Figure 4. A comparison of the Pioneer 11 FGM and VHM measurements of the 

Jovian magnetic field within 6 Rj. Quantization effects are evident at 
large radial distances, but the systematic differences between the two 
measurements are clearest near close approach. 

Figure 5. Pioneer 11 FGM inversion example at the final (third) iteration. 

Figure 6. Jovian surface isointensity contour map illustrating a 15 
eigenvector solution using the Pioneer 11 FGM observations. 

Figure 7a, b. Jovian surface isointensity contour maps obtained by perturbing 
the FGM 15 ev solution by a parameter set (proportional to the 15th ev) 
to which the observations are least sensitive. These maps illustrate the 
differences in surface field magnitude appropriate to the range of 
internal field models consistent with the observations. 

Figure 8. Pioneer 11 FGM inversion example illustrating a 13 eigenvector 
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solution. 


Figure 9. Jovian surface iso intensity contour map of the 13 ev solution using 
Pioneer 11 FGM observations. 

Figure 10. Jovian surface isointensity contour map of the 15 ev solution 

model, obtained from the Pioneer 11 FGM observations and a model magnetic 
field in which external sources are explicitly included. 

Figure 11. Jovian surface isointensity contour map of the 15 ev solution using 
Pioneer 11 VHM observations. 

Figure 12. Pioneer 11 VHM inversion example. 

Figure 13. Jovian surface isointensity contour map of a 15 ev solution using 
Pioneer 11 VHM observations and a model magnetic field in which external 
sources are explicitly included. 

Figure 14. Comparison of the GSFC 15 evs magnetic field model and the 0^ 

model at Jupiter's surface. In the upper panel the difference in field 
direction (degrees) is contoured. In the bottom panel the normalized 
difference in local field magnitude (percent) is contoured. 

Figure 15. Comparison of the JPL 15 evs magnetic field model and the SHA 23 
model, differences contoured as in Figure 14. 
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0.216 

-0.071 

0.831 


PARAfCTEHS AND PARAfCTER ASJUSTKNTS. ITERATION • 3 

HEIGHTED R«S DEVIATION 0.00530 EXTERNAL 
UEIGHTED RRS DEVIATION 0.00500 INTB9ML 


J 

9 

1 

PARAHETEK > 

4.352 

AOJ * 

0.000 

STD 

DEV 

S 

0.0072 

< 0.00211) 

LIN * 

0.7987 

J 

s 

2 

PARAfCTER » 

-0,396 

AOJ * 

0.000 

STD 

DEV 

3 

0.0060 

( 0.00287) 

LIN » 

0.7963 

J 

s 

3 

PARANETER » 

0.624 

AOJ * 

0.000 

STD 

DEV 

3 

0.0101 

( 0.00233) 

LIN = 

0.75S9 

vl 

S 

4 

PIASTER » 

-0.197 

ADJ 3 

0.000 

STD 

DEV 

3 

0.0107 

( 0.00433) 

LIN = 

0.6245 

J 

s 

5 

PARABETER » 

■^.6K 

ADJ > 

0.000 

STD 

DEV 

3 

0.0105 

( 0.00332) 

LIN * 

0.5687 

J 

3 

6 

PARANE’i'ER » 

0.295 

ADJ * 

0.000 

STD 

DEV 

S 

0.0102 

( 0.00493) 

LIN * 

0.6032 

J 

s 

7 

PARANETER » 

0.015 

AOJ * 

0.000 

STD 

DEV 

S 

0.0087 

( 0.00297) 

LIN * 

0.6921 

J 

s 

8 

PARABETER * 

-0.139 

AOJ * 

0.000 

STD 

DEV 

3 

0.0084 

( 0,00444) 

LIN * 

0.5025 

J 

s 

9 

PARABETER * 

-0.138 

ADJ * 

0.000 

STD 

DEV 

S 

0.0287 

( 0,00709) 

LIN = 

0.3350 

J 

s 

10 

PWABETER * 

-0.307 

AOJ « 

0.000 

STD 

DEV 

3 

0.0129 

( 0.00319) 

LIN * 

0.3106 

J 

s 

11 

PARABETER- » 

-0.015 

ADJ * 

0.000 

STD 

DEV 

S 

0.0156 

( 0.00622) 

LIN = 

0.3421 

J 

3 

12 

PARAMETER » 

-0.312 

ADJ * 

0.000 

STD 

DEV 

3 

0.0153 

( 0.00332) 

LIN * 

0.2904 

J 

3 

13 

PARABETER » 

0.117 

ADJ * 

0.000 

STD 

DEV 

S 

0.0191 

( 0.00552) 

LIN * 

0.3176 

J 

3 

14 

PARMCTSR » 

0.194 

ADJ * 

0.000 

STD 

DEV 

3 

0.0086 

( 0.00477) 

LIN * 

0.4552 

•J 

= 

15 

PARABETER'* 

-0.122 

ttJ * 

(1.000 

STD 

DEV 

3 

0.0205 

( 0.00948) 

LIN = 

0.3636 
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***** ITERATION NUHScR 2 ii n iii i iiiiiiii n iii n ii i iiii i ii i iiiiii 


wP. PIONcSF II AT JUPITER R < 5 W 
INTERNAL 3 EXTERNAL 0 CNDER CF PIT 
m. OBSERVATIONS 
NINSER OP EIGcM/ECTORS > IS 
PMAKETcR VECTOR : 

4.064 -0.432 0.567 -0.09S -0.934 0.368 <0.101 -0.423 -0.306 

-0.596 0.523 -0.261 -0.365 0.014 0.351 

SINGULAR VMJJES OF A 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

1 487.543 383.907 370.964 249.973 230.272 176.354 I36.28I 124.343 104.962 91.721 71.262 45.689 43.120 32.646 23.903 

V fMTRIX 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

IS 

1 

0.772 

-0.235 

0.544 

0.008 

0.128 

-0.013 

0.018 

-0.022 

0.012 

-0.031 

0.032 

-0.042 

-0.15? 

0.027 

0.019 


0.257 

0.896 

0.049 

0.152 

-0.186 

-0.162 

0.075 

-0.080 

-0.124 

0.067 

-0.067 

-0.033 

-0.064 

0.042 

-0.0.30 

3 

-0.514 

0.170 

0.775 

-0.140 

0.133 

0.096 

0.146 

-0.009 

0.085 

0.043 

0.009 

-0.143 

0.047 

0.010 

0.081 

4 

-0.089 

-0.081 

0.067 

-0.282 

-0.174 

-0.745 

-0.201 

-0.396 

0.220 

-0.220 

-0.104 

0.082 

-0.048 

0.019 

0.037 

5 

-0.041 

0.055 

-0.204 

-0.197 

0.539 

-0.236 

0.529 

0.138 

-0.170 

-0.178 

0.037 

-0.069 

-0.325 

0.124 

-0.015 

6 

-0.083 

-0.023 

0.101 

0.327 

0.332 

-0.394 

-0.425 

0.351 

-0.026 

0.493 

-0.023 

0.078 

-0.078 

0.101 

-0.199 

7 

-0.117 

-0.064 

0.004 

0.822 

0.120 

-0.052 

0.164 

-0.271 

0.284 

-0.304 

0.062 

-0.064 

-0.049 

-0.003 

0.079 

8 

-0.101 

-0.205 

0.061 

0.157 

-0.589 

-0.209 

0.179 

0.449 

-0.237 

-0.099 

-0.072 

-0.323 

-0.252 

0.209 

0.081 

9 

-0.073 

-0.039 

0.112 

0.086 

-0.080 

0.085 

0.153 

0.087 

-0.064 

-0.040 

-0.572 

0.694 

-0.231 

-O.liK 

0.021 

10 

0.015 

0.037 

0.058 

-0.013 

-0.162 

-0.084 

0.199 

0.204 

0.190 

O.Oll 

0.607 

0.578 

0.133 

0.337 

0.103 

11 

-0.084 

0.073 

0.114 

0.016 

0.028 

0.059 

-0.264 

0.137 

-0.257 

-0.612 

0.136 

0.105 

0.057 

-0.009 

-0.639 

12 

-0.055 

-0.072 

0.074 

0.101 

-0.044 

-0.189 

0.041 

-0.118 

-0.521 

0.088 

0.405 

0.089 

-0.018 

-0.660 

0.188 

13 

-0.107 

0.104 

-0.050 

-0.033 

-0.011 

0.253 

-0.391 

-0.085 

0.109 

-0.057 

0.275 

-0.002 

-0.737 

0.045 

0.173 

14 

-0.079 

-0.153 

0.048 

0.082 

-0.021 

0.089 

0.032 

-0.547 

-0.538 

0.213 

-0.016 

0.098 

-0.037 

0.544 

-0.105 

15 

-0.009 

-0.090 

-0.001 

-0.037 

-0.217 

0.003 

0.353 

-0.173 

0.238 

0.361 

0.140 

-0.062 

-0.258 

-0.212 

-0.664 


PARM1ETERS AND PARAICTER ADJUSTNENTS. ITERATION • 2 

WEIGHTED RBS DEVIATION 0.00093 EXTERNAL 
WEIGHTED RItS DEVIATION 0.00126 INTERNAL 


J 

S 

1 

PWAHETEK » 

4.064 

AOJ X 

0.000 

STD 



0.0047 

( 0.00228) 

LIN X 

0.9515 

J 

X 

2 

PWAHETEH » 

-0.452 

ADJ X 

0.000 

STD 

DEV 


0.0040 

( 0.00265) 

LIN X 

0.949? 

J 

X 

3 

PARA1ETER « 

0.567 

ADJ X 

0.000 

STD 



0.0053 

< 0.00256) 

LIN X 

0.9231 

J 

X 

4 

PARMCTER « 

-0.098 

ADJ * 

0.000 

STD 



0.0071 

( 0.00563) 

LIN X 

0.7710 

J 

X 

5 

PARIVCTER 3 

-0.934 

AOJ X 

OlOOO 

STD 

DEV 


0.0103 

( 0.00525) 

LIN X 

0.7612 

J 

X 

6 

PWAHETER > 

0.363 

AOJ X 

0.000 

STD 

DEV 


0.0109 

( 0.00599) 

LIN X 

0.700? 

J 

a 

7 

PARAMETER > 

-0.101 

ADJ X 

0.000 

STD 

DEV 


0.0069 

( 0.00446) 

LIN X 

0.8384 

J 

X 

8 

PARAMETER « 

-0.423 

AOJ X 

0.000 

STD 

DEV 


0.0123 

( 0.00528) 

LIN X 

0.6721 

J 

a 

9 

PARMSTER X 

-0.306 

AOJ X 

0.000 

STD 

DEV 


0.0193 

( 0.01139) 

LIN X 

0.3029 

J 

X 

to 

PARA«TER X 

-0.596 

ADJ X 

0.000 

STD 

DEV 


0.0192 

( 0.01208) 

LIN X 

0.4544 

J 

a 

11 

PARAMETBT « 

0.528 

AOJ X 

0.000 

STD 

DEV 

a 

0.0236 

( 0.00982) 

LIN X 

0.3583 

J 

a 

12 

PARAMETER « 

-0.261 

ADJ X 

0.000 

STD 

DEV 

X 

0.0227 

< 0.01063) 

LIN X 

0.3856 

J 

a 

13 

X 

-0.365 

ADJ X 

0.000 

STD 


X 

0.0200 

( 0.00920) 

LIN X 

0,3571 

J 

a 

14 

PARAfCTER = 

0.014 

AOJ * 

0.000 

STD 


= 

0.0186 

( 0.00326) 

LIN X 

0.4452 

J * 

15 

PARAMETER » 

0.351 

AOJ X 

0.000 

STD 


X 

0.0254 

( 0.01038) 

LIN X 

0.4638 
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